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Abstract: We discuss a method for a model independent reconstruction of the CMB 
temperature fluctuation power spectrum on small and intermediate scales that is geared 
to individual experiments. The importance of off-diagonal correlations for determining the 
shape of the power spectrum is emphasized and some examples of a reconstruction method 
are given. By using this method to "map" the power spectrum on the scales to which they 
are sensitive several experiments could be combined to map out the full power spectrum, 
and test consistency of observed features. For example we find that the GUM scan of 
the 3rd flight of the MAX experiment prefers a positive slope to the power spectrum near 
£ ~ 160, which provides weak evidence for the presence of a Doppler peak. 
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1. Introduction 

The study of anisotropies in the Cosmic Microwave Background (CMB) radiation has 
now reached the level where experiments on a wide range of scales are reporting detections 
or significant upper limits. It seems appropriate therefore to consider ways of analyzing 
the data, with a view to reconstructing the power spectrum of temperature fluctuations. 
One possibility is to give a model (e.g. CDM) of the radiation power spectrum in terms 
of a few parameters and then use many experiments to constrain the model parameters. 
In this letter we will focus on the complimentary, model independent, approach in which 
a parameterization of the power spectrum is developed on a per-experiment basis and 
individual experiments are used to constrain that part of the power spectrum to which 
they are most sensitive. Since on large scales the inversion of the sky-map into a power 
spectrum is a well studied problem we will focus here on individual experiments on smaller 
angular scales. Given that experiments now cover the full range of the power spectrum 
from degree scales up (see e.g. Bond 1993, White, Scott & Silk 1994), this in principle 
allows a reconstruction of the power spectrum. 

2. Power spectrum and window functions 

In this letter we will assume that the CMB temerature fluctuations are gaussian dis- 
tributed so that all of the information resides in the correlation matrix or power spec- 
trum. It is conventional to expand the temperature fluctuations in spherical harmonics 
T = Ylitm a tmXe.m and denote the power per mode by Cf. 

( a £m a e'm') ens = Ci Sg'gSm'm (1) 

where the angled brackets represent an average over the ensemble of temperature fluctu- 
ations. The Harrison-Zel'dovich power spectrum corresponds to constant £{£ + l)Cg or 
CJ 1 oc £(£ +1). In general one expects the Cg to be a power-law for small £ with more 
structure at higher £ (smaller angular scales) . The spectrum for CDM for example, viewed 
as £(£ + l)Cg vs ln£, is a flat line until £ ~ 100 when it rises into two "Doppler" peaks 
before falling off exponentially at large £ ~ 1000 (for a discussion of power spectra in var- 
ious models see Bond & Efstathiou 1987, Holtzman 1989, Vittorio & Silk 1992, Sugiyama 
k Gouda 1992, Dodelson & Jubas 1993, Stompor 1993, Crittenden et al. 1993, Bond et 
al. 1993). 
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When fitting model parameters to data one is interesting in computing the correlation 
matrix of the temperature fluctuations and from this a likelihood function (see e.g. Read- 
head & Lawrence 1992). The correlation matrix is a sum of two parts: the experimental 
error matrix (usually assumed diagonal: (7i<%) and a "theoretical" correlation matrix 



Cth(m, n i) = < T ( n 0^(n J )) ens = ^ + 1) C t W t (m, nj) , (2) 

71 e=i 

Here T(n) is the measured temperature assigned to direction n and Wi is the window 
function for the experiment (e.g. see Srednicki & White 1994 for a discussion of window 
functions and Bond 1993, White, Scott & Silk 1994 for plots of window functions for 
existing experiments). The method proposed below does not rely on being able to construct 
Wi{p.i, n_j) and so is more general than just Eq. (2). However for simplicity we will phrase 
our discussion as if W^(rij, n^) were known, as it is for most small-scale experiments. For 
a given C tn the likelihood function is expressed in terms of CV, = C tn (nj, n^) + cr^- as 

£oc-^exp [-^C^T,] (3) 

and provides a convenient way of using all the available experimental information to con- 
strain the model (in this case C t h). 

Usually when people discuss "the" window function for a particular experiment they 
are referring to Wi(n, n) (the window function at zero-lag) appropriate for computing the 
diagonal entries of C tn . This window function defines a "bandpass" outside of which the 
experiment has no sensitivity: the small £ cutoff is set by the amplitude of the chop and the 
high £ cutoff is controlled by the beam width (see e.g. Bond 1993). Integrating the power 
through the bandpass gives AT rms , a fact which has been exploited by several authors 
(see e.g. Bond 1993, White, Scott & Silk 1994) to estimate the amplitude of the power 
spectrum measured by the many experiments on degree scales. It has long been realized 
that separating the scales of the beam chop and width increases the area under the window 
function or the sensitivity of the experiment. 

Such an increase in the width of Wi doesn't necessarily mean less resolving power for 
Ci however. Off-diagonal entries of C t h (i-e. ^ rij) potentially contain more information 
than just the integrated power. They provide a lever arm which gives sensitivity to the 
shape of the power spectrum, as is reflected in the shape of W^n^ rij) . In figure 1 we 
show Wi{9 = 0°) and Wg(6 = 2°) (cos 6 = rvrij, with points separated parallel to the chop 
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direction) for a hypothetical experiment with 2-beam, square-wave chop of amplitude 2° 
and gaussian beam of FWHM 1°. We choose = 2° since Wg(2°) is approximately the 
derivative of Wg(0°). This derivative or"dipole-like" form provides maximum sensitivity 
to the slope of the power spectrum through the bandpass (to which C t h(0) is mostly 
insensitive) . 

The diagonal elements of Cth coming from Wg(n,n) allow an estimate of the power 
through the bandpass. From the above we see that without losing sensitivity to the am- 
plitude, the off-diagonal elements provide (model independent) information on the shape 
of the power spectrum over the same range of I. 

3. Parameterizing the power spectrum 

We wish to parameterize the Cg in a simple way that is tailored to the experiment 
under consideration. Start by defining Dg = £(£ + l)Cg and x = ln£ and extend Dg to a 
smooth function D(x) in the obvious manner. Since on general physical grounds we expect 
the power spectrum to be smooth (the neighbouring Cg are highly "correlated") we can 
expand D(x) in a Taylor series about the peak of the window function, xo, as 

D{x) = D(x Q ) [l + m(x - xo) + \m'{x - x ) 2 H ] (4) 

where to and to' are related to the derivatives of the power spectrum at xq. Due to the 
finite beam width and chopping, any experiment is sensitive to the power spectrum only 
over a limited (and usually small) range of x, defined by the window function at zero-lag. 
So for the purposes of each experiment, Eq. (4) should be a good expansion. In what 
follows we will take as given that our expansion of D(x) is not to be extrapolated outside 
of the range of sensitivity of the experiment. 

One problem with this expansion is that it does not enforce the physical condition that 
D{x) > for all x. Thus the allowed region of to, to'... is constrained in a non-trivial (but 
calculable) way. In principle one could get around this by expanding not D(x) but lnD(x) 
in a power series, but this has the technical disadvantage that the expansion is no longer 
linear in D(x). An expansion in which C t h oc C t ^ (1+toC^ +• • •) is more computationally 
efficient when performing fits to the data because , ... only have to be computed 
once. This has motivated us to stick with Eq. (4). If a Monte-Carlo analysis of the data is 
used to constrain the parameters, the advantage of linearity could be overlooked in favor 
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of direct imposition of the constraint D > 0. Eventually the data should constrain the 
parameters sufficiently that this question will not be an issue. 

4. Two examples of the method 

As an illustrative example of these ideas we have taken Eq. (4) and truncated at linear 
order (which is all the current sensitivities merit) . Converting to the conventional notation 
we assume 



where £ is the peak of Wi(0°). We computed C t h(0°, m) and C t h(2°, m) using Eq. (5) and 
the window functions of Figure 2. We find C t h(0°) (the integral through the bandpass) 
is mostly insensitive to m while C t h(2°) has good discriminating power between different 
slopes. We have, normalizing at m = 0, 



Obviously the more off-diagonal elements one has with good signal-to-noise, the more in- 
formation one has about the shape of the power spectrum through the bandpass, and 
the more one can constrain higher orders in the expansion (the "shape" of the power 
spectrum). In the above C t h(2°) ~ 0.2C t h(0°) at m = which sets the scale of sensitiv- 
ity needed to extract the full m-dependence from the signal. If noise is a problem one 
could increase C t h(6>) by measuring correlations at other separations, but these would have 
reduced sensitivity to the slope of the power spectrum. 

Clearly an experiment which observes only at widely separated points will have no 
non-vanishing off-diagonal correlations and sensitivity only to the integral of the power 
spectrum though the bandpass. Alternatively an experiment which samples too closely will 
waste time recording "redundant" information and not reduce the experimental uncertainty 
on the data points sufficiently to constrain the higher moments of the power spectrum. The 
question thus arises as to what observing strategy maximizes the constraints on the power 
spectrum, analyzed in this manner. Such a question cannot be answered without also 
considering other experimental limitations, which will clearly vary from one experiment 
to another. The question can be analyzed by the individual groups in preparing their 



e(e + i)Ct = e (e + i)c io i + min 



T 



(5) 



C t h(0°; m)/C th (0°; m = 0) = 1 - 0.2m 
C th (2°;m)/C t h(2°;m = 0) = 1 - 1.7m 



(6) 
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observing strategy. In general one gains information on the power spectrum by working 
with a strategy which (as much as possible) provides strong off-diagonal signals. 

As an application of this method to a 'realistic' scenario we have computed the like- 
lihood function in (Cg ,m) space, assuming Eq. (5), for the GUM (Meinhold et al. 1993) 
scan of the 3rd flight of the MAX experiment. The data consists of 165 points in a two- 
dimensional "bow-tie" pattern, with varying spacing and a strongly anisotropic correlation 
matrix. The fitting procedure and data used are as discussed in Srednicki et al. (1994). 
For the parameters appropriate to MAX3 Wg(n, n) peaks at ln£ — 5 (£ ~ 160) with a 
FWHM of ln£ ~ 1.4. Due to the large number of points and the wide range of correlations 
available we already find a weakly significant constraint. [The MuPeg data set (Gundersen 
et al. 1993) from the same flight consisted of relatively few, equally ("widely") spaced data 
points at constant elevation, and we find no significant constraint on m.] This analysis 
provides weak evidence that GUM prefers a positive slope over its range of sensitivity. The 
contours of the likelihood function are shown in figure 2 (we have restricted m G (—1,1) 
to enforce Cg > over the range of £ to which MAX is sensitive) . For comparison, if one 
weights by (2£+ l) 2 Wf the C £ for CDM (see Eq. (2)) with h = \ and fi B = 0.01-0.10, the 
best fits to the form of Eq. (5) are m = 0.5-0.6. The low values of the slope come from the 
fact that the MAX window function is centered near the maximum of the first "Doppler" 
peak for the models considered here. 

The analysis presented above neglects the important issue of foreground contamination, 
and the data are clearly not at the level where m is severely restricted. However it serves to 
show that the method works in principle. As sensitivities increase, a more detailed analysis 
along these lines is sure to provide information on the shape of the power spectrum. A fit 
including a quadratic component, for example, would allow one to test whether the power 
spectrum exhibits a peak in the range of £ probed (at present the likelihood function for 
MAX is very broad in m' and no conclusion can be drawn) . 

5. Conclusions 

In the future, with reduced experimental errors and more coverage, the method pre- 
sented in this letter could be used to map the power spectrum, over the range of sensitivity 
of individual experiments, in a way which provides more and more shape information as 
higher off-diagonal correlations are observed above the noise. Pasting together constraints 
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of this form from many experiments should allow one to reconstruct the full power spec- 
trum in a model independent manner. 

Acknowledgements 

I would like to thank Emory Bunn, Joanne Cohn, Douglas Scott, Joe Silk and Mark 
Srednicki for many useful conversations and comments on this and related work. This 
research was supported by the NSF and the TNRLC. 

References 

Bond, J. R., 1993, in Proceedings of the IUCAA Dedication Ceremonies, ed. T. Padman- 
abhan, John Wiley & Sons, New York, in press 

Bond, J. R. & Efstathiou, G., 1987, MNRAS, 226, 655 

Bond, J. R., Crittenden, R., Davis, R. L., Efstathiou, G., Steinhardt, P. J., 1993, preprint, 
Penn 

Crittenden, R., Davis, R. L. & Steinhardt, P. J., 1993, Ap. J., 417, L13 
Dodelson, S. & Jubas, J. M., 1993, Phys. Rev. Lett, 70, 2224 
Gundersen, J. O., et al., 1993, ApJ, 413, LI 
Holtzman, J. A., 1989, ApJ (Supp), 71, 1 
Meinhold, P. R., et al., 1993, ApJ, 409, LI 

Readhead, A. C. S. & Lawrence, C. R., 1992, Ann Rev Astron & Astrophys, 30, 653 

Srednicki, M., White, M., Scott, D. & Bunn, E., 1994, Phys Rev Lett, to appear 

Srednicki, M. & White, M., 1994, CfPA preprint, in preparation 

Stompor, R., 1993, Astron Astrophys, in press 

Sugiyama, N. & Gouda, N., 1992, Prog. Theor. Phys., 88, 803 

Vittorio, N. & Silk, J., 1992, ApJ, 385, L9 

White, M., Scott, D. & Silk, J., 1994, Ann Rev Astron & Astrophys, in press 

7 



Figure Caption 



Fig. 1. The window functions W £ (6 = 0°) (solid) and W £ (9 = 2°) (dashed) for a hypothet- 
ical experiment with 2-beam, square wave chop of amplitude 2° and a gaussian beam of 
FWHM 1°. 

Fig. 2. The contours of the likelihood function in (C £o ,m) for the MAX-GUM data. The 
contours are in units of \o. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/astro-ph/9402038vl 
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This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/astro-ph/9402038vl 



